Effect of the tetrahedral distortion on the electronic properties of iron-pnictides 
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We study the dependence of the electronic structure of iron pnictides on the angle formed by 
the arsenic-iron bonds. Within a Slater-Koster tight binding model which captures the correct 
symmetry properties of the bands, we show that the density of states and the band structure are 
sensitive to the distortion of the tetrahedral environment of the iron atoms. This sensitivity is 
extremely strong in a two-orbital (d^, d yz ) model due to the formation of a flat band around the 
Fermi level. Inclusion of the d xy orbital destroys the flat band while keeping a considerable angle 
dependence in the band structure. 

PACS numbers: 75.10.Jm, 75.10.Lp, 75.30.Ds 



The recent discovery of high temperature supercon- 
ductivity in iron arsenides 1 has triggered off intense re- 
search in the condensed matter community. The charac- 
teristic building blocks of these materials are FeAs lay- 
ers with each Fe atom surrounded by four arsenic atoms 
in a distorted tetrahedral geometry. Fe atoms form a 
square lattice, while As atoms located at the center of 
each square are displaced above and below from the Fe- 
plane in a checkerboard form (see Fig. [1]). Undoped 
compounds are compensated semimetals which undergo 
a structural and a magnetic transition^* 3 - Upon chemical 
doping, the structural distortion and magnetic order dis- 
appear and superconductivity emerges. Superconductiv- 
ity has been reported so far in three families: rare-earth 
oxyarsenides 1 ReFeAsO, also denoted 1111 systems, with 
Re a lanthanide atom, TFe2As2 with T an alkaline earth 
metal and two FeAs planes per unit cell, 4 called 122, and 
LiFeAs. 5 Iron arsenides can be considered part of a larger 
family, the iron pnictides, with As substituted by an iso- 
valent pnictogen element, as P. 6 

Hopping between Fe atoms via As is expected to give 
the dominant contribution to the kinetic energy^ and to 
the exchange interaction^ Therefore some dependence of 
the electronic properties on the angles formed by the Fe- 
As bonds, i.e on the distortion of the As tetrahedron, 
is expected. This dependence is of crucial importance 
for both weak coupling models ^ 10 ! 11 based on nesting 
properties, and strong coupling models^ based on su- 
perexchange interaction, aimed to explain the properties 
of these systems. For each family of iron-pnictides, the 
As tetrahedron presents different distortions. It is almost 
regular— in the 122 family with the four As-Fe-As angles 
close to the ideal value 109.47°, while it is elongated in 
LiFeAs, 5 with As atoms further from the Fe-plane than in 
the regular tetrahedron (the As^-Fe-As* 025 angle equals 
102.8°). Here As top refer to an As atom placed above the 
Fe plane. In the 1111 family the tetrahedron is squashed 
and the As £ °^-Fe-As^ or P to P-Fe-P top angle depends on 
specific composition, being equal to 113.7° in LaFeAsO 1 
and 120.6° in LaFePO. 6 Differences in the electronic 
properties of these compounds have been attributed to 
the different distorsion of the tetrahedra. 1 - The As-Fe- 




FIG. 1: Lattice structure of the FeAs layers with each Fe 
atom surrounded by four As atoms in a distorted tetrahedral 
geometry. The Fe atoms form a square lattice. At the center 
of each square, above and below in a checkerboard form, lie 
the As atoms, a, the angle formed by the Fe-As bond and the 
Fe-plane, is related to the As top -Fe-As top (or Fe-As-Fe, with 
the Fe atoms being next nearest neighbors) angle r/ as 2a = 
180° — rj. For instance, for a regular tetrahedron r\ — 109.47° 
and a = 35.26°. 



As angle is also sensitive to doping» i 14 i 15 i 16 i 17 and can be 
deeply modified under pressured An example of the lat- 
ter is the collapsed tetragonal phase found in CaFe2As2 
under pressured In Refs. Il3l andll4l a possible correlation 
between the superconducting critical temperature and 
the As-Fe-As angle has been suggested with deviations 
within a family from the regular tetrahedron being detri- 
mental for superconductivity. In Bai_ ;E K ;E Fe2As2 at op- 
timal doping the tetrahedron changes from squashed (in 
underdoped) to elongated (in overdoped).— It has been 
argue d 13 i 14 that there is a correlation between stronger 
interaction and higher critical temperature in iron pnic- 
tides due to a narrower bandwidth for the regular tetra- 
hedron. This argument seems to be supported by the 
lower critical temperatures and lack of structural distor- 
tion and magnetism in LiFeAs and LaFePO. An unusu- 
ally large sensitivity of the iron moment 19 and the band 
structur o 1Q i 2Q i 21 to the separation of the As atoms with 
respect to the plane has been also found in Density Func- 
tional theory calculations. 

In this paper we analyze the angle dependence of the 
band structure within a tight-binding model in which 
hopping between Fe atoms is assumed to be mediated 
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by As and its magnitude is calculated within the Slater 
Koster framework. 22 We show that the two-orbital model 
which only includes d xz and d yz is extremely sensitive to 
changes in the angle formed by the As-Fe bonds. Such 
a sensitivity originates in the appearance of flat bands 
in the (0,0) - (±tt, 0) and (0,0) - (0, ±tt) directions for 
the regular tetrahedron case. The corresponding peak in 
the density of states is found at the Fermi level at half 
filling. The flat bands disappear when the tetrahedron 
is distorted. Within the present model this result is in- 
dependent of any fitting parameter. Inclusion of the d xy 
orbital modifies this picture but maintains the depen- 
dence of the density of states and band structure on the 
angle. 

Two- orbital model. We construct a tight binding model 
to describe the band structure of the FeAs layers around 
the Fermi level in the tetragonal phase. We describe 
the distortion of the tetrahedron in terms of the angle a 
formed by the Fe-As bond direction with the Fe-plane, 
see Fig. [TJ The advantage of using this angle, instead of 
the more common Fe- As-Fe or As- Fe-As ones, is that it is 
uniquely defined. On the contrary, there are two different 
Fe- As-Fe angles (depending on whether the Fe atoms are 
nearest or next-nearest neighbors) and, for non-regular 
tetrahedra there are also two different As-Fe- As angles 
(between top As and between top and down As). The 
correspondence between a and the Fe- As-Fe or As- Fe-As 
angles is straightforward (see figure caption of Fig. [1]). 

According to first principles calculations several bands 
cross the Fermi level and the density of states (DOS) 
at the Fermi level is dominated by Fe-d orbitals with 
little weight of pnictogen p-orbitals j 7 i 2Q i 23 We work in 
the Fe-square lattice and take x and y axes along the Fe- 
Fe bonds. Only the Fe-orbitals are included in the tight 
binding. As atoms only enter in the model indirectly via 
the Fe-Fe hopping amplitudes. Under these assumptions 
the Hamiltonian is given by: 

+ ^/3,7 C i,j,/3,cr C «+ 1 5 i+ 1 5 7,c r + tptf C i,j,P,a C i+lj-l,'y,o' 



+ h.C.)+ ^ e P C h,f3,* C ij,(3,<r ~ V 



(i) 



Here z, j and a label the sites in the 2D lattice and the 
spin respectively, while /?, 7 refer to the Fe-d orbitals in- 
cluded in the tight binding, e is the on-site energy, and 
\i is the chemical potential. We first restrict ourselves to 
two degenerate d xz and d yz orbitals, as in other minimum 
models recently discussed j 24 i 25 Due to the degeneracy of 
d xz and d yz orbitals the onsite energies e xz = e yz only 
shift the bottom of the bands and for simplicity we take 
them equal to zero. The hopping amplitudes t are cal- 
culated from orbital overlap integrals within the Slater- 
Koster framework, 22 which captures the correct symme- 
try properties of the energy bands. Direct Fe-Fe hopping 
is neglected, i.e. we assume that all hopping takes place 
via As atoms. Neglecting the energy splitting between As 
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FIG. 2: Angle dependence of the Fe-Fe hopping parameters 
(left) and Fe-As overlap integrals (right) which enter in the 
two-orbital model. Fe-As overlap integrals and Fe-Fe hopping 
parameters are in units of pd a and pd^/(ed — e p ) respectively. 
pdj, = —0.2. Shaded areas correspond to the experimentally 
relevant angles. 



p XyV and p z orbitals and to second order in perturbation 
theory: 

%z,xz = 2 [fe,-) 2 - fe,-)' - fe,-)'] /fe - £ P ) ( 2 ) 
t V xz,xz = 2 [fe,-)' - fe,-)' + fe,-)'] /fe - *p) (3) 

and tf,z,yz = tV x Zi xzi tV yz, V z = Kz.xz' Hopping to near- 
est neighbors between xz and yz vanishes. Next near- 
est neighbor hopping amplitudes satisfy t" z = —t' xz , yz , 

t" — f' f" — f f and t f = +' with 

xz,xz xz^xzi yz,yz yz,yz aj±lu - b yz,yz xz,xz VV1U±1 

t'xZ,XZ = [{tx,XZ? + (ty, XZ ) 2 ~ (t z , xz f] /(€ d ~ € p ) (4) 



^xz,yz — [ tx,xztx,yz ty,xzty,yz tz,xz^z,yz] / {^d £p) 

(5) 

with e p and the on-site energy of the pnictogen-p and 
Fe-d orbital, and ts,p, with 5 = x,y,x and (3 = xz,yz, 
the overlap between As-p^ and Fe-dp orbitals. These ex- 
pressions take into account that two As atoms mediate 
the hopping between nearest neighbors, while only one 
As atom is involved in the hopping between next near- 
est neighbors. Overlap between As-p orbitals and Fe-d 
orbitals is given in terms of two integrals, taken as dis- 
posable constants, pd a and pd^, where as usual a and ir 
refer to the component of angular momentum around the 
axis^ The sign of pd n is taken opposite to that of pd a . 26 
The angle-dependence is evident in the Fe-As overlaps: 

t x ,xz = sin a ( ^ cos 2 a(pd a ) + sin 2 a(pd n ) ] (6) 



sin a cos' 



: a(V3(pd a )-2(pd 7r )) (7) 
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FIG. 3: Top (bottom) figures: Density of states (band struc- 
ture) for the two-orbital model corresponding to the squashed 
(left), regular (middle) and elongated (right) tetrahedron. 
The corresponding angles are a = 32.71°, 35.26° and 39.13°, 
respectively. In the top figures the green (light) lines show the 
total density of states. Black and red lines give the density 
of states corresponding to the E~ and E + bands (Eq. [9} re- 
spectively. The density of states is smoothed by a lorentzian 
of width 0.02. Energies are in units of pd^/(ed — e p ). 



tz.xz — 



COS OL 

7T 



(^/S sin 2 a(pd a ) + (1 — 2 sin 2 a)(pd^) 



(8) 

and t x ,yz — ty^xzi ty,yz — tx,xz-> and t Zt)XZ — tz,yz- The 

angle dependence of the hopping parameters which enter 
in Eq. (pQ) and of the Fe-As overlaps are shown in Fig. O 
A strong dependence is seen, specially in the Fe-Fe hop- 
ping, for the angles found experimentally (shaded area in 
Fig. [2]). In this figure Fe-Fe hopping parameters and Fe- 
As overlaps are given in units of (pd a ) 2 / (e p — e^) and pd G 
respectively. In these units all the hopping amplitudes 
and overlap integrals depend on a unique free parameter 
pd n . 

Diagonalization of the Hamiltonian results in two 
bands 



E ± (k) = c+(k) - (i ± Je_(k)2 + eUk) (9) 



with 



e +( k ) = (^L,^+^L,*J( cos k x +cos k y )+4t' z cos k x cos k y 

(10) 

e -( k ) = (txz,xz - ty xz,xz)( COsk x - COS ky) (11) 



e xy (k) = 4t f yZiXZ sin k x sin k y 

The density of states (DOS) for a squashed (a = 32.71°), 
a regular (a = 35.26°) and an elongated tetrahedron 
(a = 39.13°), corresponding to As^-Fe-As^ angles 
of 114.58°, 109.47° and 101.74° respectively, and pd^ = 
—0.2 is plotted in Fig. [3] (top panels). In this figure the 
chemical potential is at co = 0, assumming half- filling. 
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FIG. 4: Energy contour plots for the squashed (top), regular 
(middle) and elongated (bottom) tetrahedral geometries for 
the same parameters as in Fig. [3] Left figures correspond to 
E~ and right figures to E + , see Eq. [9] The Fermi surface 
are in black (darker lines). The topology of the Fermi surface 
changes with the angle. 



Two peaks are observed in the DOS for the three an- 
gles. The DOS of each of the two bands E ± (Eq. [9]) is 
also plotted in these figures. It can be seen that each of 
the peaks comes from one of the two bands. The height 
of the high energy peak is very strongly sensitive to a. 
In particular, the intensity of this peak is much higher 
for a regular tetrahedron. Inspection of the band struc- 
ture in the lower plots of Fig. [3] reveals the existence of 
a fiat band along (0, 0) — (7r, 0) in this case. By symme- 
try, the band is flat also along the (0,0) — (0, ±tt) and 
(0,0) — (— 7r,0) directions. The flat bands in the regular 
tetrahedron geometry (a = 35.26°) appear for any value 
of pda and pd^. From the expressions for the hopping 
amplitudes given above, for this a it can be shown that 
txz,xz + 2t' xz xz = 0, which cancels the k x dependence of 
the band structure for k y = (alternatively, it cancels 
the /c y -dispersion for k x = 0). This equality for the reg- 
ular tetrahedron comes from t XjXZ = t ZyXZ , see Eqs. (|6j) 
and ([8]) and right panel in Fig. 2. The topology of the 
energy contours and, therefore, that of the Fermi surface, 
is also strongly sensitive to a as shown in Fig. [H for the 
same parameters used in Fig. [3j 

Due to the change in the energy contour (and the Fermi 
surface) topology, and the very high peak in the DOS, if 
interactions were included, the magnetic and supercon- 
ducting properties of this two-orbital model would be 
extremely sensitive to the angle a. This is particularly 
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FIG. 5: Density of states for the three orbital model corre- 
sponding to the squashed, regular and elongated for the same 
angles as in Fig. [3] pd^ = —0.2 and e xy — e xz , yz = 0.15. The 
density of states is smoothed by a gaussian of width 0.03. 
Energies are in units of pd 2 a / (ed — e p ). 

evident in a picture which ascribes the magnetic transi- 
tion to a spin density wave due to nesting. We show now 
that including a third orbital maintains some angle de- 
pendence in the electronic properties although when the 
hybridization with the third orbital is strong the condi- 
tions for the flat band are not fulfilled. 

Three orbital model According to LDA£ i 2Q i 23 the d^ 
orbital contributes to the Fermi surface and it is ex- 
pected to be relevant in any description of iron pnic- 
tides. 27 Hamiltonian (pQ) can be straightforwardly gen- 
eralized to include the d xy or any of the e g orbitals*^ 
In the case of the regular tetrahedron, including only the 
crystal field due to the arsenic environment, the three t2 g 
orbitals (d xy , d yz , and d xz ) are degenerated The inclu- 
sion of the crystal field produced by the Fe environment 
and/or the distortion of the tetrahedra breaks this de- 
generacy. As the As atom approaches the Fe-plane the 
d xy orbital becomes higher in energy as compared to the 
(still degenerate) d xz and d yz orbitals. The level splitting 
£xy — £xz,yz introduces a new parameter in the calculation. 
When both Fe and As crystal field effects are included 
£xy — £xz,yz is slightly larger for the squashed tetrahe- 
dron than for the regular or elongated one. Expansion of 
hopping parameters to second order in perturbation the- 
ory also contributes to the renormalization of the on-site 
energies e xy and e xz , yz . 

For large level splittings e xy — e xz , yz > W. with W the 
bandwidth of the two-orbital model, hybridization be- 



tween d xy and the bands E + and E~ of the two-orbital 
model is weak and the picture described above remains. 
However, in real materials the level splitting is expected 
to be smaller. 20 Fig. [5] shows the DOS of the three or- 
bital model for the same geometries as in Fig. [3] and with 
£xy — £xz,yz = 0.15, and pd^ = —0.2. For simplicity, we 
neglect the angle-dependence of the level splitting. Fig. [5] 
reveals that a clear dependence of the DOS on the angle 
a still prevails within the three-orbital model. Note that 
the Fermi level is expected to be around the center of 
the band where the effect of changing the angle is very 
strong. 28 However, the origin of the dependence of the 
DOS on the angle a is qualitatively different to the one 
found in the two-orbital model: Due to the hybridiza- 
tion with the d xy orbital the flat band that appeared in 
the regular tetrahedron geometry within the two-orbital 
model is partially destroyed in the three-band model. 

In conclusion, we have analyzed the angle dependence 
of the band structure within a Slater-Koster based tight 
binding model for iron pnictides. In both two and three 
orbital models the density of states and therefore the 
Fermi surface are strongly sensitive to changes in the 
angle formed between the Fe-As bonds. This result is 
relevant for both weak coupling theories based on the 
nesting mechanism, which is very sensitive to the shape 
of the Fermi surface, and for strong coupling theories, 
since superexchange also depends on orbital overlaps. In 
the two-orbital model a flat band is formed when the As- 
tetrahedron is regular what results in a strong peak in the 
density of states. This result is robust against changes 
in the fitting parameters. Hybridization with a third or- 
bital modifies this picture although the dependence on 
the angle remains, in particular for low energies. In our 
description, the bandwidth does not seem to depend very 
strongly on the angle, though it seems to be slightly nar- 
rower for a more elongated tetrahedron. 
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